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ABSTRACT 

Procrustes rotation involves fitting a factor 
pattern matrix to a specified target matrix in factor analysis. These 
rotations ats useful for the investigator who wishes \o see how well 
his data can be made to fit a hypothesised factor pattern matrix. The 
mathematical problems involved in these transfoi nations are outlined 
and computer algorithms for both oblique and orthogonal solutions are 
presented. (C K > 
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Abstract 

Procrustes rotation involves fitting a factor 
pattern matrix to a specified target matrix. This 
report briefly outlines the mathematical problems 
involved and provides computer algorithms for both 
oblique and orthogonal solutions. 
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Those familiar with Creek mythology will remember that Procrustes was 
a character whom Theeeus encountered in his extended travels. Procrustes 
took pride in having beds which fitted all travelers: those persons who were 
tco tall he cut down to size and those who were too short he stretched. The 
name, Procrustes, is particularly appropriate for the factorial rotation 
methods to be discussed here. 

It sometimes happens in the course of research involving factor analysis 
thai* an investigator raa., want to see how well his data can be made to fit 
a hypothesized factor pattern matrix. The problem may be stated more 
formally as: Given a factor pattern matrix. A, and a hypothesized factor 

pattern matrix, B, how may a transformation matrix, T, be obtained such that 

AT = B + E 

is satisfied with the restrictions that 

tr (E'E) = minimum 

and 
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Horst (1956) seems to have been one of the first to propose a workable 
solution for this. Ahmavaara (1957) considered the problem, as did Hurley 
and Cattell (1J02). All of these authors gave essentially the same 
oblique solution. 

The solution is very easily computed. Solving 

T*= (A ’A) _1 A'B 

and normalizing the columns of T*gives T, the desired trar.sf ornat ion matrix. 
Multiplying to obtain AT gives a matrix which is a least squares fit of 
B with the restriction that 

Vi = 

The above solution is oblique. An orthogonal solution utilizes the 
additional constraint that 

T f T « TT’ = I. 

Thi3 somewhat more difficult problem was examined in derail by Sch*6nemann 
(I960). His method begins by taking the minor product of A and B as 

S = A’B. 

then calculating S f S and S3'. Obtaining all roots and vectors of these two 
matrices, preferably by a Jacobi method, gives 

S'S a VDV 1 

and 

3S' - WTO 1 . 

A probiem occurs at thin stage concerning the proper orientation of the 
orthogonal vectors in V and W. This is quite easily solved by computing 

W'SV = 
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and checking ?ach diagonal element in . If the diagonal element is nega- 
tive, the corresponding column of W is reflected. 

With the new reflected matrix, W*, we may proceed to obtain 

T = W*V' 

as the desired orthogonal transformation matrix. Calculating AT gives an 
orthogonal least squares fit of B, 
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PROCRUC 



Program Setup Instructions 



A. Heading Card -- Any alphanumeric message to be printed at be- 
ginning of output * 



B. Problem Card 
Column 
i - 7 
8 -12 

11-14 

15-16 



Item 

Punch "PROBLEM". 

Iteration tolerance for Jacaim 
Root and vector Subroutine. 

(Try 00100 or 00010.) 

Number of variables in each matiix 

Number of factors in each matrix 
(Note — the loading matrix and the 
target matrix must be the same size) 

C. Input Variable Format Card. -- card with any F-type format. 

Both input matrices must have the same format. 

D. Date~-Must be on cards. Matrices must have variables in rov^s, 
factors in columns. Input two complete matrices in the following 
order. 

1. Factor loading matrix to be transformed. 

2. Target matrix. 

This is a repeating program. Any number of problems may be run by 
repeating steps A through D» 

E. Finish — At the end of all problems place a card with "FINISH 11 
punched in columns 1-6. 
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RUNX(S) . 

LGO. • 

-00 

PROGRAM PI<0CI<US< INPUT .OUTPUT. TAPE 5=I*VUT.TAPE6=0UTPUT ) 
DIMENSION A <25.1^1 «B 1 25 » 10 » • AT < 25 » 10 ) » T ( 10 • 10 J 
*.HDG(8) * F N". ' t 8 ) 

C THIS PROGRAM WAS WRITTEN BY CARL JtNSEMA IN JANUARY . 1971. 

C IT PROVIDES BOTH OtiLIOUE AND ORTHOGONAL PROCRUSTES SOLUTIONS 

C THE METHODS COMt FROM PS YCHOMET R I KA < MARCH. 1 966 . P . 1 -1 0 ) . 

C A= INPUT LOADING MATRIX. 

C B<= TARGET MATRIX, 

C T = TRANSFORMATION MATRIX. 

C AT* 3 TRANSFORMED LOADINGS. 

C TOL= ITERATION TOLERANCE. 

C NV= NUMBER OF VARIABLES. 

C NF = NUMBER OF FACTORS. 

10 READ < 5 • 1 1 HOG 1 PROM 1 TOL » NV . NF »FMT 

1 FORMAT < 8A10/A7 » F 5 • 5. 21 2/8A 1 0 ) 

I F I PROB. NE.7HPRUBLEM) GO TO 10000 
IF (HDG 1 1 t »E0. 10HFINISH I GO TO 10000 
I F (PROD . EO . 7HFINISH ) GO TO 10000 
I F < F MT C 1 ) .EO.IOHFINISH ) GO TO 10000 
WR I TE ( 6 « 2 I HDG. TOL .NV.NF.FMT 

2 FORMAT! A5H1UNIVERSITY OF WASHINGTON BUREAU OF TESTING 

1//SAH PROCRUSTES TRANSFORMATION PROGRAM// 1 X . 8 A 10 
?///2«H iteration TOLERANCE 10X.F9.5 

1// 2 OH NUMBER OF VARIABLES 10X.I3 

4//20H NUMBER OF FACTORS 10X.I3 

S///15H INPUT FORMAT 8 A 10 I 

WRI TE< 6. 3 I 

3 FORMAT ( 5 I / ) » 1 5H INPUT LOADINGS ) 

DO 20 I * 1 1 NV 

READ (6 *FMT I ( A ( I . J ) . J= 1 »NF > 

WRITE! 6» A) I*IAU* Jl.JaltNF) 

A FORMAT!// AH ROW 1 3 » ( /1 5F9 . 3 > > 

20 CONTINUE 

WR I TF 1 6 . 5 I 

5 FORMAT ( 5 ( / ) 1 1AH TARGET MATRIX) 

DO 3o 1 = 1. NV 

READ IS .FMT HOI I . J) . J = 1 .NF ) 

WRI TE(6> A I I » t P. C I » J) >0*1 »NFI 

30 CONTINUE 

DO SO 1=1 »2 

IF(I.EO.l) CALL UBLIPROI A.B.AT.T.TOL.NV.NF) 

IFU.EQ.21 CALL OR ThPRO ( A . B . AT . T , TOL .NV . NF ) 

WRITE (6. 6) 

6 FORMAT(5(/).22H TRANSFORMATION MATRIX) 

00 AO J= 1 .NF 

WRI TE(6»A)J.!Tt J.K) .X=1.NF) 

AO CONTINUE 

WRI TEI6.7I 

7 . • F0RMAT!S(/).21H TRANSFORMED LOADINGS ) 

DO J= 1 .NV 

WRI TEI6»A I J« <AT ( J. X) .X=l »NF ) 

50 CONTINUE 

GO TO 10 
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10000 WRl TEI 6 ♦ 8 ) 

8 FORMAT ( 10 l / ) » 1?H END OF J03. ) 
STOP 
END 



SU3 ROUT I NE UBL I FRO ( A , B » AT , T » T OL, N V * MF ) 

C OBLIQUE PROCRUSTES 

DIMENSION A ( 25 , 10) .13125 * 10) .AA( 20.10 )»AB( 10.10) ,T( JO, lO) > 

♦AT ( 25.10) »E ( 10} 

WR I TE » 6 * 1 ) 

1 FOWMAT(lOI/)»lX»2l(]H*l//20H OBLlOUE SOLUTION // ,lX , 2 1 ( 1H* J > 

C OBTAIN MINOR PRODUCT MOMENT OF LOADING MATRIX. 

00 10 I = 1 .NF 
DO i 0 J=1 .NF 
AA ( I . J ) = 0 .0 
DO 1 0 K=1 »NV 

AA< I » JliAAII » J)+A( <» I I*A(K. J) 

10 continue 

C OBTAIN ROOTS AND VECTORS. 

CALL JACSIM(AA»E.TOL?NF1 

r ROOTS ARE IN DIAGONAL OF UPPER PART OF AA AND ALSO IN E. 

C EIGNVECTOKS ARE COLUMNS IN LOWER PART OF AA. 

C INVERT MATRIX. 

DO 2C 1=1. NF 
E( I 1=1 .O/E ( I > 

20 CONTINUE 

DO 30 1=1 .NF 
I NF = I +NF 
DO 30 J= 1 »NF 
JNF= J+NF 
AA ( I i J ) =0 . 0 
DO 30 K= ) .NF 

AA( I . J) = AAII . J)+AAI INF»K)*E(IA)*AA( JNF.O 
30 CONTINUE 

C MULTIPLY THE LOADING MATRIX BY THE TARGET MATRIX 

DO AO I = ] .NF 
DO H J=1«NF 
AB ( I . J)=0.0 
DO Ao K = 1 . N V 
A8i I .J)=A3t I .JJ+AIK.I 
AO CONTINUE 

C MULTIPLY AA BY At) TO GET TRANSFORM MATRIX. 

. DO 5 0 1 = 1 iNF 
DO 50 Je i .NF 
T ( I ,J)=0.0 
00 50 K=1 »NF 
T( I .J)=T( I .JJ+AAM 
50 CONTINUE 

C NORMALIZE THE TRANSFORMATION MATRIX BY COLUMNS 
DO 70 J = 1 .NF 
SUMSOaO.O 
DO 60 1 =1 .NF 
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SUMsO=$UMSQ+ ( T < I i J } **? ) 

60 CONTINUE 

SUMSQ=SQRT(SUM50) 

DO 70 l = 1 »NF 
7 i I » J 1 *T ( I , J) /SUMSQ 
70 CONTINUE 

C MULTIPLY LOADINGS BY TRANSFORM 

00 80 I = 1 »NV 
DO 80 J=1»NF 
AT ( I »J)=0.0 
DO 80 K= 1 »NF 

AT( I »J J=AT( I ,J)+A( I »KI«TCK. J) 
80 CONTINUE 

RETURN 
END 



SUBROUTINE OKTHPRO ( A ♦ B. A T . T » TOL ,N V » NF ) 

C ORTHOGONAL PROCRUSTES * SCHOiNEMANN METHOD! PSYCHOMETKI<A»MARCH, 1966 ) 
DIMENSION A( 25 . 10) »R ! 25 . 10 ) ♦ B AAR l 20 . 10 ) »ARBA(20»10 ) » 

*AR(i O, Iiy > »T( 10, 10) rAT I25»}0( »E(10 ) 

WRITEI6.1) 

l FORMAT ( 1 0 ( / I *1X*21 OH*)//? OH ORTHOGONAL SOLUTION //IX ,21 ( JH« > > 

C OBTAIN MINOR PRODUCT OF LOADINGS AND TARGET. 

DO 10 I =1 »NF 
DO i 0 J=1 .NF 
ABt I ♦ J ) =0 . 0 
DO 1 0 K=l iNV 

ARI I » J ) = AB ( I , J J + A( K* J )*BIX*J) 

10 CONTINUE 

C OBTAIN MINOR AND MAJOR PRODUCTS OF A6. 

DO 20 1 = 1. NF 
DO 20 J=1 i NF 
BAAdU >J)=0.0 
ABBA t ! » J ) = C . 0 
DO 20 K=I»NF 

BAAd ( I , J)=BAAR( I , J > + AB ( K , I > * AB < K. J ) 

ABBA! I ♦ J ) = ABBA.I I , Jj+ABII »K)*AB<U»K) 

20 CONTINUE 

C OBTAIN ROOTS AND VFCTORS 

CALL JACSI MI OAAB ,E ♦ TOL *NF ) 

CALL JACSI M (ABBA.E > TOL »NF ) 

C VECTORS A«E IN LOWER PART OF THE MATRICES. 

C BEGIN REFLECTION OF ABBA VECTORS. 

C (SEE P • 8 > BOTTOM PARAGRAPH OF SCHONEMANN ARTICLE) 

C MULTIPLY TO OBTAIN WSV. 

DO 30 1 = 1, NF 
DO 30 J = 1 »NF 
ABBA ( I , J) =0.0 
DO 30 <*l*NF 
XNF=K+NF 

ABBA! I » J ) = ABBA I I , J ) +ABRA I KNF , 1 )*ARU»J) 

30 CONTINUE 

o 
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DO 4 0 1 = 1* NF 
DO 40 J=] »NF 
AB ( I » J ) =0 . 0 
DO 4 0 K = 1 » N," 
knf=k+nf 

AB< I ♦ J)=A8I f t J ) + A3BA ( 1 > < ) *BAAB l KN F ♦ J \ 

40 CONTINUE 

C CHEC< 0 1 AGONAL OF W SV» REFLECT COLUMN IF DIAGONAL ELEMENT IS NEGATIVE. 

DO 5C I =1 »,NF 

IF! ABU . IT.GE.O. ) GO TO 50 

DO SO J=1 * NF 

JNF=J+NF 

ABBA ( JNF ♦ I ) = Aflf)A ( JNF » I > * l “1 * > 

50 CONTINUE 

C END OF REFLECTION 

C OBTAIN TRANSFORMATION MATRIX BY MULTIPLYING VECTORS. 

DO 60 I =1 »NF 

inf=i+nf 
DO 60 J=1 »NF 
JNF=J+NF 
T ( 1 1 J) *0.0 
DO 60 K = 1 »NF 

T ( I . J > =T ( l . JJ+ABBA ( INF *K>*BAAN(JNF*K) 

60 CONTINUE 

c multiply loadings qy transform. 

DO 70 l * 1 i NV 
DO 70 j«l »NF 
AT ( I *J>=0.0 
00 70 K*1 »NF 

AT! I *J)=AT! I • J)+A! I ,KI*T(K.J) 

70 CONTINUE 

RETURN 
END 

SUBROUTINE JACS I M ( R , D » P * N > 

c this subroutinf calculates all roots ano vectors cf a matrix 

C USING A JACOBI METHOD. 

DIMENSION R! 20.10) .0! 1 0 ) 

N1=N+1 
N1 1--N-1 
N2=N *2 

DO 1U I =N 1 » NO 
DO i 0 J= 1 »N 
R( I . J) =0.0 
10 CONTINUE 

00 20 1=1, N 
N I *N + I 
R(NI.!)=1.0 
20 CONTINUE 

DO 90 L=) *100 
00 NO I = 1 » N 
D! I ) =R ( I ♦ I ) 

DO CONTINUE 

DO 5 0 I =1 ,N) 1 

n-i+i 

DO 50 J»IltN 
DR = RU *1 )-R! J,J) 
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A-SQRT(DR**?+4.*R( I » J>**2> 
A=SQRT ( < A + DR »/ (2 ,*A) ) 

8=SQRT ( 1 .-A**2 ) 

C = S I GN < 1 . *R( I iJl) 

DO 4 0 K=? >N? 

U=R (Kill *A*C+R f 6 t J) *B 
R(<,J)=~R<K.»I )*0*C + R(K, J>*A 
R ( K » I )=U 
40 CONTINUE 

DO 50 K= 1 » N 
U=R ( I , K) *A*C+R { J *K ) *B 
R(J »<)=-«< I ,<l#tJ*C+R{J*K|*A 
Rt l »K. I =U 
50 CONTINUE 

DO 60 I =1 »N 

D< I ) =ABS { D < I)-R< I > I ) > 

60 CONTINUE 

5=0.0 

DO 70 I = 1 » N 
S=NAXl (SiDlMI 
70 CONTINUE 

DO 80 I = 1 » N 
Dt I )=R t I , II 
80 CONTINUE 

IF('-.-P) 100,100,90 
90 CONTINUE 

IOy return 

END 
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